Create map of CA counties

Mutate Census Population data for 2020.

%%Mutate Census Population data for 2010. %%```{r} dec_vars_2010 <- listCensusMetadata( name = “2010/dec/pl”, type = “variables” )

smc_pop_2010 <- getCensus( name = “dec/pl”, vintage = 2010, region = "block:*“, regionin =”state:06+county:081“, vars =”group(P1)" ) %>% mutate( block = paste0(state,county,tract,block), pop = P1_001N )


Isolate the blocks for North Fair Oaks 2010.
%%``` {r}
library(devtools)
library(tigris)
smc_blocks_2010 <- blocks("CA", "San Mateo", year = 2010, progress_bar = F)
nfo_boundary <- places("CA", progress_bar = F) %>% 
  filter(NAME == "North Fair Oaks")
nfo_pop_2010 <- smc_pop_2010 %>% 
  left_join(smc_blocks_2010 %>% select(block = GEOID10)) %>% 
  st_as_sf() %>% 
  st_centroid() %>% 
  .[nfo_boundary, ] %>% 
  st_set_geometry(NULL) %>% 
  left_join(smc_blocks_2010 %>% select(block = GEOID10)) %>% 
  st_as_sf()

Isolate the blocks for North Fair Oaks 2020.

Place these isolated blocks into a map.

Overlap between 2010 and 2020

%%{r} leaflet() %>% addProviderTiles(providers$CartoDB.Positron) %>% addPolygons( data = nfo_pop_2020, stroke = F, fillOpacity = 0.5 ) %>% addPolygons( data = nfo_pop_2010, color = "red", weight = 0.75, fill = F )